Phase transitions in the steady state behavior of mechanically perturbed spin glasses 
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We analyze the steady state regime of systems interpolating between spin glasses and ferromag- 
nets under a tapping dynamics recently introduced by analogy with the dynamics of mechanically 
perturbed granular media. A crossover from a second order to first order ferromagnetic transition as 
a function of the spin coupling distribution is found. The flat measure over blocked states introduced 
by Edwards for granular media is used to explain this scenario. Annealed calculations of the Ed- 
wards entropy are shown to qualitatively explain the nature of the phase transitions. A Monte-Carlo 
construction of the Edwards measure confirms that this explanation is also quantitatively accurate. 

PACS: 05.20 y Classical statistical mechanics, 75.10 Nr Spin glasses and other random models, 
81.05 Rm Porous materials; granular materials 
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In complex systems such as granular media the en- 
ergy available due to thermal fluctuations is not suffi- 
cient to cause particle rearrangement, hence in the ab- 
sence of external perturbations the system is trapped in 
a metastable state. A granular media may be shaken me- 
chanically and experiments reveal a steady state regime 
defined by an asymptotic density |^ . The non trivial be- 
havior of these systems, such as slow relaxation dynamics 
and hysteresis effects arises from the fact that such sys- 
tems have an extensive entropy of metastable or blocked 
states. A vertically tapped system of hard spheres tends 
to a random close packing |^], whereas a horizontally 
shaken system crystallizes |^], the stationary states ob- 
tained in these systems are not theoretically understood 
for the moment. Edwards has proposed that one 
may construct a thermodynamics over metastable states 
in the same way as Boltzmann and Gibbs developed over 
microstates, his hypothesis is that the equilibrium mea- 
sure is flat over all blocked or metastable states satisfying 
the relevant macroscopic constraints. This is an ergodic 
hypothesis which is conceivably true in the case of ex- 
tensive manipulations such as stirring, pouring, shaking 
an tapping. This scenario had been recently tested on a 
wide range of models by comparing dynamical time aver- 
aged quantities in the steady state with those predicted 
by a flat measure which is either calculated numerically 
or analytically . In a number of systems the Edwards 
measure was seen to work extremely well, although, re- 
assuringly given the strength of the hypothesis, there are 
cases where it does not appear to be applicable [|j . 

Given the recent interest in the Edwards measure and 
the fundamental and industrial importance of mechani- 
cally driven systems it is natural to examine phase tran- 
sitions in this context. As pointed out by Edwards and 
coworkers, the construction of a steady state measure 
for mechanically perturbed systems would allow one to 
predict, for example, whether two powders or granular 



media form a mixture or phase separate upon stirring or 
shaking. In mean field implementations of the Edwards 
measure to binary granular systems second order phase 
transitions have been found between the mixed and sep- 
arated phases . 

Recently Q there has been interest in using spin sys- 
tems to examine tapping like dynamics on systems with 
an extensive entropy of metastable states. These sys- 
tems, though far from being realistic models of granular 
media, provide a testing ground for the development of 
thermodynamic like theories for granular systems and in- 
deed seem to show much of the phenomenology seen in 
real systems. 

The models we shall consider are Ising spin systems 
on random thin graphs. A random thin graph is a col- 
lection of N points, each point being linked to exactly 
c of its neighbors, c therefore being the connectivity of 
the graph. The spin glass/ferromagnet model we shall 
consider has the Hamiltonian H — 
where the Si are Ising spins, is equal to one if the 
sites i and j are connected and zero otherwise. The 
fact that the local connectivity is fixed as c imposes 
the local constraints Uij ~ c, for all sites i. For 
the following, we have fixed c — 3. The couplings 
Jij are independent random variables with distribution: 
P{Jij) = a5{Jij + 1) + (1 — a)S{Jij — 1) (a being a pa- 
rameter in [0, 1/2]). 

We define the tapping dynamics as follows: In between 
the taps, the system has a natural zero temperature re- 
laxational dynamics, it evolves under a random sequen- 
tial single spin flip dynamics where only moves which re- 
duce the energy are allowed. When blocked, it is tapped 
with strength p G [0,1/2], that is to say each spin is 
flipped with a probability p, the updating at this point 
being parallel - this corresponds to the extensive manip- 
ulation. The system is then evolved by the zero tem- 
perature dynamics until it becomes once again stuck and 
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the tapping is repeated. After a transient regime, the 
system reaches a stationary regime characterized by an 
energy per spin E(p). Numerical simulations show that 
E{p) is a monotonically decreasing function of p - the 
lighter one taps, the lower the energy obtained. Here we 
make the analogy with experiments on vertically tapped 
granular media , where (in the reversible part of the ex- 
periment) the asymptotic compactivity increases as the 
system is more lightly tapped. Our numerical simula- 
tions are reversible and one can move up and down the 
curve of E{p) if the system is tapped for a sufficiently 
long time after changing p. Here the macroscopic quan- 
tity which is fixed on average by the tapping dynamics 
is E{p), whereas in granular material it is the average 
volume V = V/N per particle. 

The canonical form of the Edwards hypothesis [0J|] 
then leads to the consideration of the partition function 
Z(/3) = J dEdm exp{-N{l3E - SEdwiE,m))), where 
SEdw{E, m) the Edwards entropy per spin at average en- 
ergy E and magnetization m per spin. Here [3 — [3{jp) is 
a Lagrange multiplier fixing E, numerical tapping simu- 
lations thus suggest that /3(p) decreases as p is increased. 
In spin glasses the Edwards entropy per spin at fixed av- 
erage energy can be calculated (or approximated) and 
can be shown to be extensive {i.e. SEdw{E) is of order 1 
and independent of N in the thermodynamic limit). As 
mentioned previously granular systems are expected to 
share this property. Let us emphasize at this point that 
this rather ambitious goal of developing a thermodynam- 
ics for granular systems does use additional information 
about the dynamics as compared to usual thermodynam- 
ics. The definition of the entropy contains the informa- 
tion about mctastability. This is a partial, though rather 
limited, use of dynamics to describe statics. 

If we assume that I3(p) is continuous on reducing 
p one expects the possibility of an onset of ferro- 
magnetic ordering, which would correspond to phase 
separation in a granular material. The Edwards 
free energy per spin /(/3) = miriE,mf{E,m, (3) = 
wnEf{E,(3) with J{E,m, (3) = E ~ SEdwiE,m)/ P and 
f{E, (3) = f{E, m{E), P), where m{E) is the point where 
f{E, m, P{p)) is a minimum with respect to m. The value 
of m{E) is thus given by dsEdw{E , m{E)) / dm = 0. In 
the models studied here we find that there exists a value 
Ec such that for E > Ec one has m{E) = and for 
E < Ec the solution m — Q becomes a local maximum of 
f{E, m, (3) and there is a second order phase transition 
with m becoming non zero. 

For E < Ec 'va the neighborhood of Ec we have found 
that, depending on the parameters of the disorder, there 
are two possible behaviors: 

(1) d'^SEdw{E,m{E))/dE'^\^^^- < in which case the 
ferromagnetic phase is stable and the system exhibits a 
second order phase transition in the energy with a dis- 
continuity in d(3/dE. 



(2) d^SEdw{E,m{E))/dE^\j^^^- > and the ferro- 
magnetic phase is unstable. In this case there is a second 
minimum of f{E) (case (2a)) or end point of f{E) (where 
the Edwards entropy vanishes, that is to say the ground 
state Eos) (case (2b)) at an energy E* which is strictly 
lower than Ec. Hence there is a first order phase transi- 
tion with a non zero spontaneous magnetization m{E*). 
We find in this case that the paramagnetic phase may 
be in fact only metastable for a range of values (3 < (3c 
i.e. the second minimum exists below Pc and may have 
a lower Edwards free energy. Comparing our simulations 
with the annealed calculation we find it is only when 
this metastable phase becomes unstable (at the spinodal 
point) that the transition occurs. This suggests that the 
tapping dynamics is not efficient at tunneling over free 
energy barriers - it is interesting to note that the fluffy 
(arch rich) phase in granular media exhibit this strong 
mctastability 

With respect to the relaxational part of the dynamics 
introduced above, the total number of metastable states 
for a given set of Jij and is given by : 

Nms{E, m) = Tr [] I ^ .h^n^jS^Sj | (1) 
1=1 j 

6{H - NE) 5{Y^ - Nm) 

i 

where 9{x) = if a; < and 6{x) = 1 if x > 0. 
The above formula expresses the fact that each spin is 
aligned with its local molecular field in a metastable 
state and hence none can flip under the relaxational dy- 
namics. The Edwards entropy of the metastable states 
of energy per spin E and at magnetization m is thus 
SEdw{E,m) = jj {{In NMs{E,m))), the brackets denot- 
ing the average over the disorder. Unfortunately this av- 
erage is not easily taken analytically at low energy, and 
we have to compute an annealed average, which will give 
an upper bound for the Edwards entropy: 

SEdUE,m)= lim ^\n{{NMs{E,m))) (2) 

In the the SK model and for Ising spin systems 
on random thin graphs H, there exists an energy Eann 
above which this annealed approximation is exact. 

The computation of this quantity is a generalization of 
that for a = 1^. One finds SEdwiE,m) is given by the 
saddle point value of the action 

S{u, V, t, P; E, m) = -/3 (f - £;) + f ln(l - 2a) (3) 
- f In [(1 - a){2t + u'^ + fv'^) - a{2tuv + 1 -|- i^)] 
+ i±™ In (/(ue'5)) -t- In f{veP)) 
_ i±™ln(l±™)-i^ln(i^) 

with respect to the parameters u, w, t and /?. The mag- 
netization of the metastable states at energy E is given 
by the maximization of (ra) with respect to m: 
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In order to compute the Edwards entropy, we solve the 
stationarity conditions coming from the extremization 
with respect to u,v,t,P of S{u,v,t, P; E,'m{E)), with 
m{E) replaced by the expression (^. The resulting value 
for P is P = ^£^^LE)^ which is the definition of the in- 
verse Edwards temperature. 

The saddle point can be found numerically to compute 
f{E, P). In the annealed approximation when a < ac — 
0.06 we find that we are in scenario (2), for a > etc we 
are m scenario (1). Shown in Figs. (|l]), (|) and (|) is the 
behavior of f{E, P) for a = 0, a = 0.04 and a = 0.12 re- 
spectively and one sees that one has the behaviors (2b), 
(2a) and (1). In the range Po < P < Pc shown in Figs. 
(|l|) and (|^) the higher energy, paramagnetic, minimum is 
only metastable. 

To test the above predictions, we have tapped systems 
with a = 0.002,0.005,0.02,0.05 and recorded both the 
energy E(p) and the magnetization m(p) in the quasi 
equilibrium regime. Systems were of size A'^ = 10^ and 
we have checked that this was large enough to consider a 
single sample instead of several realizations of disorder. 
The curves for E{p) are shown in Fig. (||). For a — 0.005 
and a = 0.002, there is a first order phase transition, 
whereas for a = 0.05 and a — 0.02 this transition is a 
second order one, in qualitative agreement with the pre- 
vious discussion of the role of the Edwards entropy. Of 
course, as the annealed entropy is not expected to be ex- 
act a low energy, the values of ac and Ef,{a) which come 
from this calculation are only approximate. Indeed, it is 
found in the tapping simulations that 0.005 < ac < 0.02, 
and E{pc{cty^ , a) differs systematically with Ec{a) found 
fromEq. @. 

In order to test this scenario in a quantitative fash- 
ion, we have computed the Edwards entropy as was done 
in One introduces an auxiliary temperature 1/Paux, 
which is a Lagrange multiplier to fix the number of spins 
which are not in the same direction as their local field. 
Then, one fixes P and does a Monte-Carlo simulation 
using the Metropolis algorithm where a randomly cho- 
sen spin is flipped with probability min (l, e~'^""^^'^«"=") , 

with PauxT-iaux = PU + PauxY^ i^^ {^iY. j^i'^ij -^ij ^ ■ 

Starting from Paux = 0, we increased it very slowly un- 
til the system was blocked in a metastable state. Re- 
peated several times, this allowed us to sample E{P) for 
systems of = 10^. At the end of the simulation, we 
have recorded both energy and magnetization. The com- 
parisons between m{E{p)) from tapping simulations and 
m{E{P)) from Monte-Carlo simulations for a — 0.05 are 
shown in Fig.(^, the agreement is excellent. In the high 
temperature region P G [0,/3c(0)], for all a E [0,1/2], 
the value of E{P) measured in the simulations were in 
excellent agreement with the results of the annealed cal- 
culation, indicating the latter is correct in this range. 



For Pc{0) < P < Pc{a) (here we mean Pc from the an- 
nealed approximation) the annealed calculation is close 
but certainly deviates from the numerically measured one 
(see Fig (^). Moreover, the value of Ec{a) found in this 
Monte-Carlo procedure was much closer to E{pc{a)~^) 
than that found in the annealed approximation. This 
means that this approximation is not exact for E when 
m{E) ^ and explains the quantitative disagreement 
between E{pc{a)~ ,a) and E*{a) obtained from (||). 

For a < a^ the simulated annealing again predicts 
Ec{a) very close to E{pc{a)^ ,a) determined by the tap- 
ping simulations. The behavior of the low energy magne- 
tization for E < E* (a) is also the same for the two sim- 
ulations, although the region of energy close to E*{a) is 
difficult to sample with the simulated annealing. Fig.(0) 
is the histogram of the energy obtained during the tap- 
ping dynamics in the stationary regime for a — 0.002 
for different p near pc, for a relatively small system size. 
There are clearly two peaks separated by a gap and the 
spread towards negative energy of the one of highest en- 
ergy becomes large when p approaches Pc, thus indicating 
the approach of the spinodal point. It is interesting to 
compare Fig. (^ qualitatively with Fig.(0), the compar- 
ison is quite striking. For P slightly smaller than Pc we 
see that the form of the high energy peak, with a small 
spread to the right and a large one to the left is predicted 
by the geometry of the first minimum of f{E,P) in Fig. 
(^ . In addition for /? slightly smaller than Pc we see that 
the ferromagnetic minimum is close to its spinodal point, 
thus explaining the large, two sided in this case, spread 
in the histogram in this region. 

In conclusion we have examined a dynamics driven by 
an extensive external perturbation on a family of sys- 
tems having an extensive Edwards entropy of blocked 
states. The existence of first and second order ferro- 
magnetic transitions are qualitatively explained by invok- 
ing the Edwards measure calculated within an annealed 
approximation. The quantitative predictions using this 
measure in Monte-Carlo simulations are also in excellent 
agreement with tapping simulations on the system. We 
have found similarly promising results in models of higher 
connectivity and in the infinitely connected SK spin glass 
|]l2| . This work suggests the possibility of using the Ed- 
wards measure to explain phase transitions in realistic 
granular systems, the main obstacle is the technical dif- 
ficulty in the computation of the Edwards measure. 
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Figure captions 

Fig. 1. Free energy f{E,l3) versus E for different vahies of 13 with a = 0. Between A and B, the local minimum 
has the lowest free energy, between B and C, it is only metastable and above C it is no longer a minimum. 

Fig. 2. Free energy f{E,f3) versus E for different values of (3 and a = 0.04 < ac. One sees that when (3 increases 
the minimum of high energy starting as a global minimum becomes metastable and then disappears at the transition. 

Fig. 3. Free energy f{E,f3) versus E for different values of (3 and a = 0.12 > ac- Here there is only one local 
minimum which moves continuously with f3. The vertical dashed lines indicates the energy E where f{E,f3) is non 
analytic. 

Fig. 4. Energy per spin E{p) in the steady state regime versus p for a = 0.002, 0.005, 0.02, 0.05, 0.5 and N = 10^ 
spins. 

Fig. 5. Comparison of the measurements of m{E) in the tapping simulations and Monte-Carlo simulations. In the 
former (a) m{p) is plotted versus E{p) and in the latter (b) m(/3) is plotted versus E{(3). 

Fig. 6. Comparison of E{P) for a = 0.05 obtained from the annealed calculation (a) and the Monte-Carlo simulation 
(b). At the right of the vertical dashed line, the annealed calculation is not exact any more. It predicts a first order 
phase transition, whereas the Monte Carlo predicts a second order one, with Ec in very good agreement with that of 
the tapping at p = 0.05. 

Fig. 7. Histogram of the energy during the tapping simulations in the steady state regime for a = 0.002 and 
N = 25000 spins during 10^ taps and for p = 0.2584, 0.2585, 0.2586, 0.2587, 0.2588, 0.2589 (reading left to right from 
the top). 
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